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Abstract 

We study ground-state phase diagrams and excitation spectra of Bose-Bose mixtures in an optical 
lattice by applying the Gutzwiller approximation to the two-component Bose-Hubbard model. A 
case of equal hoppings and equal intra-component interactions for both components is considered. 
Due to the existence of inter-component interaction, there appear several quantum phases, such as 
the superfluid, paired superfluid, and counterflow superfluid phases. We find that the transition 
from superfluid (SF) to Mott insulator (MI) with even filling factors can be of the first order for 
a wide range of the chemical potential. We calculate the excitation spectra as a useful probe to 
identify the quantum phases and the SF-to-MI transitions. In the excitation spectra of the SF 
phase, there are two gapless modes and a few gapful modes, which respectively correspond to 
phase and amplitude fluctuations of the order parameters. At the SF-to-MI transition point, we 
show that the energy gaps of certain amplitude modes reach zero for the second-order transition 
while they remain finite for the first-order one. Since the excitation spectrum can be measured 
by the Bragg scattering, we calculate the dynamical structure factor by using the linear response 
theory. We consider three types of density fluctuations, and show that the density fluctuations are 
coupled to different excitation branches in different quantum phases. 

PACS numbers: 67.60.Bc, 03.75.Lm, 05.30.Jp, 05.30.Rt 
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I. INTRODUCTION 



Since the observation of the quantum phase transition from the superfluid (SF) to Mott 
insulator (MI), ultracold atoms trapped in an optical lattice have provided unique fields for 
studies of strongly correlated quantum matter l|,[2|. In particular, the experimental creation 
of mixtures of two types of bosons in optical lattices j^-^j has opened up new possibilities 
to explore exotic quantum phases. Previous studies have analyzed the two-species Bose- 
Hubbard model, which quantitatively describes systems of Bose-Bose mixtures in optical 
lattices, and predicted various quantum phases, including SF, MI, phase separation, paired 
superfluid (PSF), counterflow superfluid (CFSF), density wave, and supersolid phases js)- 
341 ]. The PSF phase is a superfluid of composite particles that consist of two bosons o 
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is present for attractive inter-component interaction 
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32] . The CFSF is expected to emerge inside the Mott insulator 



for repulsive inter-component interaction, and is interpreted as a superfluid of composite 
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es consisting of a boson of one type and a hole of the other 
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321 ] . In the PSF (CFSF) phase, only the in-phase (out-of-phase) motional degrees 



of freedom exhibits superfluidity, i.e. the motion is dissipationless, while the out-of-phase 
(in-phase) motion is prohibited. Despite the intensive interest, these two exotic phases have 
not been experimentally observed so far. 

Hu et al. have reported the dipole oscillations in the SF, PSF, and CFSF phases in an 
one-dimensional optical lattice using the time-evolving block decimation method 32j- In 
order to understand such a dynamical behavior, it is important to reveal the details of the 
excitations of Bose-Bose mixtures in an optical lattice. The excitations of a one-component 
Bose-Hubbard model 
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36] have been well understood from previous theoretical and ex- 



perimental studies. In the MI phase, the lowest two branches of the excitation spectrum are 



gapful and correspond to the particle- and hole-excitation modes 



37H40L In the SF phase, 



the excitation spectrum has one gapless mode and gapful modes 39 
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451 ] . The gapless 



mode corresponds to oscillations of the phase of the superfluid order parameter and is known 
as the Bogoliubov mode while the lowest gapful mode is regarded as the amplitude mode 
in the vicinity of the SF-to-MI transition at commensurate fillings. It is also well-known 



that the amplitude mode becomes gapless at the critical point. 



been experimentally detected via the Bragg spectroscopy [46l-l52| and the lattice-amplitude 



'he excitation spectra have 
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modulation 
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55| and been successfully used to characterize the two phases and the phase 



transitions. Likewise, it is expected that elementary excitations can be used to distinguish 
the different phases and the SF-to-MI transitions in Bose-Bose mixtures. 

In this paper, we study the quantum phase transitions and the excitation properties 
of Bose-Bose mixtures at zero temperature in an optical lattice using the Gutzwiller ap- 
proximation (GA). We assume equal hoppings and equal intra-component interactions for 
both components. To describe the SF and MI phases, we directly apply the GA to the 
two-component Bose-Hubbard model. However, this approach fails to capture the PSF and 
CFSF phases because these phases arise from the hopping of particle pairs or anti-pairs. 
In order to analyze the phases, we use the effective Hamiltonian describing the degrees of 
freedom of pairs or anti-pairs within the second order-perturbation theory. First, we deter- 
mine the ground-state phase diagrams of this system, for which our main focus is placed 
on the first-order SF-to-MI transition. Kuklov et al. have performed Monte Carlo simula- 
tions on so-called two-component J-current model, which is a classical counterpart of the 
two-component Bose-Hubbard model at even total fillings, and found that the SF-to-MI 
phase transition at even total fillings can be of the first order for a certain region of the 
inter-component interaction [l4|. More recently, Chen et al. have applied a method based 
on the tensor product states to the two-component Bose-Hubbard model with attractive 
inter-component interaction and the hardcore constraint, and shown that the transition 



from incommensurate SF to MI can be of the first order as well [30|] . In the present paper, 
we will show that this is also the case for the two-component softcore Bose-Hubbard model 
with repulsive inter-component interaction and that the region of the first-order transition is 
extended to a wide range of the chemical potential. Second, in order to study the excitations, 



we extend the method for calculating excitation spectra of the Bose-Hubbard model [431. 156 1 
to the system of Bose-Bose mixtures. We discuss the effect of inter-component interaction 
on the excitation spectra in the SF and MI phases. We show that the gap of the in-phase 
amplitude mode vanishes at the SF-to-MI transition for odd total filling and another am- 
plitude mode also becomes gapless at the second-order SF-to-MI transition for even total 
filling. Finally, we investigate the response to density perturbation by applying the linear 
response theory to the lattice system. In the SF, PSF, and CFSF phases, we calculate the 
dynamical structure factors regarding responses to the three types of density perturbation, 
namely one-component, in-phase, and out-of-phase fluctuations, in order to show that these 
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phases can be identified by means of the Bragg scattering techniques. 

This paper is organized as follows: in Sec. Ull we explain formulations, which consist of 
the two-component Bose-Hubbard model, GA, effective Hamiltonian, linearized equation of 
motion and the linear response theory. These formulations are used in Sec. IHIllTVl and El In 
Sec. IIIH we obtain the phase diagrams for several parameters. We discuss first-order phase 
transitions from SF to MI with even total fillings for repulsive inter-component interaction. 
In Sec. HV] we determine excitation spectra in the MI, SF, PSF, and CFSF phases. We see 
several changes in the excitation spectra when hopping amplitude decreases from SF to MI. 
In Sec. El we apply the linear response theory to a Bose-Bose mixture system, and discuss 
the response to the Bragg scattering. The conclusions are presented in Sec. ED 



II. FORMULATION 



A. Two-component Bose-Hubbard model 



We consider a system of a D-dimensional hypercubic optical lattice loaded with a mixture 



of two different types of bosons, which can be two different hyperfine states 



7 



8[, atomic 



species {4], and isotopes [57[. Assuming that the optical lattice is sufficiently deep com- 
pared to the chemical potential, we model the system by the two-component Bose-Hubbard 
model [361 ] . 



a=l,2 (jj) % 



n n _i[n 
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where b ai (b a ,i) is the creation (annihilation) operator of the a(= 1 or 2) component at site 
i, n a ^ = b ai b at i is the number operator, and (i, j) denotes the sum over the nearest neighbor 
sites. t a , U a , and fi a are the hopping amplitude, the on-site intra-component interaction 
and the chemical potential, respectively; U12 is the on-site inter-component interaction. In 

58 
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experiments the magnitude of U12 can be controlled by the Feshbach resonance 
and by the component-dependent optical lattice §[ 0] . We determine the quantum phases 



by the compressibility k = n t 1 (f^J , the polarizability <p = (j^j > an d the order 
parameters <3> Qii = (b a ,i), $f = (&i,t&2,i), and $| = where n t 
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finite 
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PSF 
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finite 





CFSF 





+ 








finite 


No LRO 


















TABLE I: List of the quantities characterizing the different phases, namely the superfiuid (SF), the 
paired superfiuid (PSF), and the counterflow superfiuid (CFSF), and the phase with no long-range 
order (LRO). + or finite means that the corresponding quantity takes the positive or finite value. 
We drop the site index i of the order parameters by assuming the spatial homogeneity. 



A*t — A*i + fi>2, A* a = A*i — fi>2, and n a is the number of particles per site (filling factor) of the 
component a. The conditions for identifying each phase are summarized in Table [H Notice 
that the incompressible phases, namely CFSF and one with no long-range order (LRO), are 
regarded as Mott insulators. Instabilities toward the phase separation and the collapse are 
characterized by the conditions <fi < and k < 0, respectively. In this paper, we consider 
the special case of U\ — U2 — U > 0, t± — £2 = t, ^1 — ^2 — A*, and \Un\ < U. If U12 does 
not satisfy this condition, the mixtures lead to the phase separation for repulsive interaction 



or the collapse for attractive interaction 



Ml, 



61] 



B. Gutzwiller approximation 

We use the GA to investigate the ground states and the excitation spectra of Eq. ([T]) at 
zero temperature. We assume the Gutzwiller-type variational wave-function written as 

00 

l*> = II E fS,n 2 (r)\n u n2h, (2) 

l 711,712=0 

where 712)1 is the Fock state with ri\ and ri2 particles of the components a = 1 and 2 at site 
i, and the variational factor fnln 2 has to satisfy the normalization condition ^2 ni „ 2 \fni,n 2 \ 2 = 
1. One can obtain the equation of motion for fn}, n2 by imposing the stationary condition 
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on the effective action J dr(^f\ihj- — H\*&), which leads to 



ih 
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II 2 



dr 



U 



U a {n a — 1) — pLUa 
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where the superfiuid order parameters for each component are $ij = 

En 1 ,n 2 /ni-i,n 2 v /7I i/S«2 and $ 2j = £ ni ,„ a /iSXa-i V^f®** ■ We define a s sum 

over neighboring sites of site i. In the ground state, the wavefunction is stationary so that 
the coefficients can be written by 



fW (t) = f w e 

Jni,7i2\ / J ni,n2 



(4) 



where fn},n 2 is the coefficient of stationary state that does not depend on time, and the 



phase factor Ui is given by 

oo 



ni,n 2 



— ni(ni - 1) - Ami + —n 2 (n 2 - 1) - jtm 2 + U\ 2 n x n 2 



J ni,n 2 



(5) 



where $ M = E m , n2 /m-i^V^/Sk and $ 2 ,i = E ni ,„ 2 /£,^-iV^/"i!«a are the superfiuid 
order parameters for stationary state. We calculate the ground state coefficients by solving 
the Gutzwiller equation with the imaginary time propagation method 62|, |63|. More 
specifically, we follow the steps as shown below; (i) take the imaginary time r' = ir and set 
adequate initial coefficients, (ii) calculate the order parameter <3>j, (iii) put $j into Eq. (J3J), 
calculate new coefficients, (iv) iterate the steps (ii) and (iii) until the coefficients and the 
average energy converge. In the actual calculation, the Hamiltonian matrix is truncated at 
finite values n c % and n c2 . In the present work, we set the initial coefficient as a real number, 
and n c \ = n c2 = n c and use a sufficiently large n c so that the results do not depend on n c . 



C. Effective Hamiltonian 



One cannot describe the PSF and CFSF phases with the equation of motion (3), because 
there the hopping term is treated as the first-order perturbation while the hopping of pairs 
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(anti-pairs), which are the essential degrees of freedom in the PSF (CFSF) phase, is a process 
of the second order with respect to the hopping. For instance, if one tries to use Eq. (3), 
these phases are not present in any regions of the phase diagrams. Moreover, it fails to 
resolve dispersion of the gapless mode in the excitation spectra. Hence, to describe the 
CFSF and PSF phases, we use an effective Hamiltonian that is restricted in the low-energy 
subspace of pairs or anti-pairs and account for tunneling of pairs or anti-pairs within the 
second-order perturbation theory. 

In the limit of t/U <C 1 , the Hamiltonian of Eq. ([1]) can be written as H = Hq + tV by 
treating the hopping term as a perturbation, where the non-perturbative Hamiltonian H 
and perturbation V are given as 



V 



E 



2 Tl>a,i\W>a,i 



E [«LA 



- 1) - fin a . 
h.c. 



(6) 
(?) 



Using t 
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re second order perturbation theory, we can derive the effective Hamiltonian 



19. 
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31] . In this paper, for simplicity, we consider the case that the amplitude of the 



inter- component is close to intra-component interaction, |C/ 12 | < U . 

In the case of attractive inter-component interaction [7 12 < 0, the low-energy subspace of 
pairs is described by a product over single-site Fock states with equal occupation of the two 
species. For < (U + U\2)jU <C 1, the effective Hamiltonian is given by 14 . 120I \l 



H, 



off 



H 2t V 



hi{hj + 1) + hj(hi + 1) + b\,p 2i bijb 2 j + fcy^ii^i 



it it 



(8) 



The appropriate Gutzwiller wave-function for the effective Hamiltonian Eq. (jHJ) is given by 



i*p>=nE/« (i) Ki 



(9) 



This wavefunction leads to the equations of motion, 

m jL /n P« = [ Un (n _ 1) _ 2 ^n + U 12 n 2 ] ft® 



2t 2 



U 
2t 2 



(10) 
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where = Yin, [fnf-i n ifnf > ]'i an d Uj = is the average particle number, which 

corresponds to the number of pairs in this phase. In the ground state, we can describe the 
coefficients /n^n 2 (r) = fn^e~ luJ ^ T , where uf is given by 



hwf = ^2 { Un ( n - 1) - 2 ^ n + u - 



12" 



2^ 



^[n(n j + l)+ri j (n + l)] 



<i)i 



f" 



p(0 



<i>» 



(11) 



where $ P = J^ n . {fnl-* n ifni^) is the pair superfluid order parameter for the stationary 
state. 

On the other hand, for the repulsive inter-component interaction C7i 2 > 0, the low-energy 
subspace of particle-hole pairs is described over single-site Fock states with uniform total 
on-site occupation n t . For < (U — U\2)/U <C 1, the effective Hamiltonian is given by 



u 



^ ^ [n a ,i{n a j + 1) + h a ,j{n a ,i + 1)] - jr ^ (^L^i^Ai + &ii&2j&2t&iij 



and the appropriate Gutzwiller wave function is 



(12) 



(13) 



From the calculations similar to the attractive case, we obtain the Gutzwiller equation, 
'^fn {l) = {I [2n 2 + nt (n t - 2n - 1)] - /xn t + C/ 12 n(n t - n)| /°« 



2^ 

77 



[2n% + n t (n t - n - + 1)] /«« 



" ^ E [^V"("t-" + l)/S + <&r>/(n + l)(nt-n)/, 



c(i) 
ra+1 



(14) 



where = J^ n . /° _* y / nj (n t — nj + l)/n^ • The ground state is described as /, 



c(i) 



F c(i) 
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fn^e 1UJ ^ T , where the phase factor u\ is given by 



— {2n 2 + rit (n t - 2n - 1)} - /m t + U 12 n(n t - n) 



2t 2 



— jj £ {2nn j + n t (n t - n - rij + 1)} 



(j>« 



c(i) 



9+2 , 



(15) 



Here $^ = 2~2 nj \f^-i'\/ n j( n t ~ n i + tyfnfj * s ^ ne counterflow superfluid order parameter 
for the stationary state. 



D. Linearized equations of motion 

In this section, we derive the linearized equations of motion, which allow us to calculate 
the energies and the wave functions of elementary excitations. In Sec. IIV[ we will use the 
derived equations to calculate the excitation spectra for several different phases. 

We consider a small fluctuation from the stationary state given by 

(16) 

Assuming that the stationary state is homogeneous, we expand the fluctuation in terms of 
the plane wave, 

SfS,n 2 (r) = £ (u nun2 ^-^) _ v* nun2M e-^-^) , (17) 

k 

where rj is the position vector of site %. Substituting Eq. (fT7|) into Eq. fl3]) and linearizing 
the equations with respect to the small fluctuations, we obtain 




where Uk and Vk are (n c + l) 2 -dimensional vectors with the components w nii n 2 ,k and v ni , n2 ^ 



f(i) I 
J ni,ri2 v 



fW +5f w (t) 
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respectively. The matrix elements of and are given as, 



.4 



(rci,"2),(n'i,n2) 



u 



-n n (n. 



1) - jin a + U\2n\n 2 - ftoJi 



&ni ,17,^ &ri2,r 



Zt f^lV^l^ni-l.n'i + $ *VM- Itfni+l.ni 



(ni,n,2),(n' 1 ,ri2) 



rt2 — ljJij 

+ $2V^2 + lOna+l,' 

- e(k) [C^ + Ci w + C 2 x + C^.J , 
e(k) \p\_-y + Di^ + Dl _ x + D!^] , 



(19) 



(20) 



r 1 

u l,m 

r 2 

n 1 



D 2 



'n' x + (1 + /)/2 v / n 1 + (1 + m)/2f* n , + , ./ Bl+ 



m,n 2 ; 



Vn' 2 + (1 + 0/2^2 + (1 + m)/2/ n * ;>n , + J n 

i,»j+m) 

i^n; + (1 + 0/2^ + (1 + m)/2f n[+w J n 
Vn 2 + (1 + /)/2v/n 2 + (1 + m)/2f n> , 2+ J nu 



(21) 



where we define e(k) = 2£ cos(fc^a), a is the lattice spacing, and D is the spatial dimen- 
sion. We can calculate the excitation spectrum of a given stationary state by diagonalizing 
Eq. (HHP. 

Since Eq. ( Tl8l) does not capture correctly the low-lying excitations of the PSF and CFSF 
phases, we also derive the linearized equations of motion from the effective Hamiltonians 
Eqs. (IH]) and ( 1T2|) using Eqs. ( ITUi) and ( fl4l) . respectively. The linearized equations of motion 
for the effective Hamiltonians Eqs. flH]) and f fT2|) take the same form as Eq. ffT8]) . For the 
former case, the basis vectors Uk and Vk are (n c +l) -dimensional vectors with the components 
u n k and v n k, and the matrix elements are given by 



A n,n 



Un(n - 1) - 2fm + U 12 n 2 



2zt 2 



2zt 2 



[n (n + 1) + n (n + 1)] - M)f > 5 nin , 



U 



$ p * (n + 1) <W,„> + $ p n^_! , 
- e p (k) [(n' + 1) (n + 1) + n'nflUfU 

eP(k) [n' (n + 1) /^/^ + (n' + 1) nj*_i#H , 

c/ 2^z= 

and Vk are (n t + l)-dimensional vectors with the components u n ^ and i> n k, and the matrix 



-r>n.n 



(22) 
(23) 



where we have defined e p (k) = jf^2 l=1 cos (kia). For the latter case, the basis vectors Uk 
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elements are given by 



a n,n 



nn.n 



|^ [2n 2 + n t (m -2n- 1)] - /m t + U 12 n(n t - 



n 



2zt 2 
2zt 2 



U 



[2nhj + n t (n t — n — hj + 2)] — hu^ j S n: n' 
5> c V( n + !) ( n t ~ «)<Wi,n' + $V n K - n + l)<J„_i,„/ 



- e c (k) \y/{n> + 1) (nt - W) (n + 1) (n* - n)f c J + J c n+1 
+ y/n~' (n t - n' + 1) n (n t - n + , 



= e c (k) [y^ (n t - n' + 1) (n + 1) K " 



+i 



+ y/(ri + 1) (n t - ri) n(n t -n + l)Jl,_J, 



c 

ra+1 



where we define e c (k) = ^ Y^iLi cos (^z a )- 



(24) 



(25) 



E. Linear response theory 

In this section, we briefly present the linear response theory applied to the equations of 
motion, namely Eqs. fl3]), ([ID]) , and (fl4|) . which describes Bose-Bose mixtures in an optical 
lattice. More detailed derivations are shown in Appendix A. In Sec. |V] we will use the 
linear response theory to calculate the dynamical structure factors for three types of density 
fluctuations. 

We consider the time-dependent perturbation of the following form: 

H pcrt (r) = -J2 (xAe'^e^ + \*G\e^e^ , (26) 

i 

where Aj is the strength of external field that will be taken to be sufficiently small, rj is a 
small constant describing the adiabatic switching on of the perturbation at t — > — oo, and 
Gi is the local operator given by 

G{ = G 

mi , ni2, ni ,772 

|mi,m 2 )i(ni,n 2 |i. (27) 

771 ,772, mi, 7772 

The external field can be expanded in terms of plane waves as Aj = ^ k Ak ]W e *. 
The fluctuation of the expectation value of an local operator of physical interest JFJ = 

Eru.na.mi.ma F ™i,m 2 ,n u n 2 K, ™ 2 )i(™l, 7l 2 |i is given by 



= E ^,a(k,^)e-^ + X ^,Gt(k, 



ikr. 



(28) 
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Within the GA, the response function ^pg(k,w) is given by (for a detailed derivation, see 
Appendix [A} 



(0\F\u)(u\G\0) _ {0\G\v){v\F\G) 
u + ir] - <jj^ v u + ir] + uj^ v 



(29) 



where |0) is the ground state, and \v) is the v th excited state. The matrix element of the 
response function is defined as 



,n 2 / j 



(0|O|^) — ^/ni,n 2 ^™i.«2,ni i «2 ti ni,n2,k V m ,n 2 ,\fl n i .™2,"i ,«2 /«i 

(z/|O|0) = ^u nin2k O niin2jniin2 f nijri2 — f ni> n 2 O nijn2>ni>n2 v^ n2 jA , (30) 



ni,n. 2 



where u 



ni,n 2 ,k 



and t> 



ni,n2,k 



are the solutions of the linearized equation of motion f[T8]) . In the 



calculations of Sec. El we set rj/U = 10~ 2 for the SF phase and r]/U = 10~ 4 in the PSF and 
CFSF phases. We choose the operators F = G and calculate the response function using 
the results of Sec. IIVBI We obtain the dynamical structure factor Sp(k,u) from following 



relation 



64| 



S p (k,u)) = Im (xpffcu)/'*) . 



(31) 



Note that this formulation is valid only at zero temperature and Sp(k, u) vanishes for u < 0. 
We normalize the dynamical structure factor as SV(k, uj) = Sir(k, oj)/Sf(^) by using the 
static structure factor given by 



S>(k) 



duSp(k, oj) 



(32) 



III. GROUND STATES AND FIRST ORDER TRANSITIONS 



In this section, we use the GA described in the previous section to obtain the phase 
diagrams in the (zt/U, fi/U) plane for equal hoppings, equal intra-component interactions, 
and several values of Uu/U, where z is the coordination number. Since we will investigate 
properties of elementary excitations in different quantum phases in Sees. IIVI and |V] it is 
useful to locate these phases beforehand. In previous studies, similar phase diagrams have 
been calculated by using mean-field theories [ll|, Q and the strong-coupling expansion 



techniques 



29| | . However, the Mott insulators with odd total fillings have not been addressed 



13 



in these studies. Moreover, while the first-order SF-to-MI transition was predicted to occur 

n 

at even total fillings [14|, it has not been confirmed that the first-order property remains for 
the transition between MI and incommensurate SF except in the case of hardcore bosons 



with attractive inter-component interaction 



30| . Below we will address these points. 



A. Phase diagrams 



In Figs. QJa) and (b), we depict the phase diagrams for repulsive inter-component inter- 
action. We see that there is a region of an incompressible MI for each integer value of the 
total filling factor. The existence of the MI phases with odd fillings has been predicted in 

nnn 

Refs. [9|, [12J, [13(. In the atomic limit (t/U = 0), the size of the MI regions in the //-axis for 
odd and even total fillings is U12 and U, respectively, and the odd-filling MI regions vanish 
at U12/U — > [17]. Once the regions of the MI phases are located, we solve Eq. (TBI) in 
order to check whether or not the CFSF order is present. It is shown that the CFSF order 
emerges in the MI phases with odd total fillings while it is not present in MI with even 
total fillings. This result is consistent with that in Refs. 

ACT 
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25 



26|. Within the 



mean-field approximation used here, the local state of the ground state at odd total filling is 

75 Ofr 1 ' fr 1 ) + fr 1 )) while that at even total fillin g is If' f >■ As indicated by the 

dashed line in Fig. QJa), when U12 is close to U, there is a wide region of the first-order phase 
transition to MI with n t = 2 on the phase boundary. We note that similar first-order tran- 
sitions have been found in multi-component Bose-Hubbard systems with inter-component 



exchange interactions 



65- 



In Sec. IIIIBt we discuss the first-order transitions in detail. 
In Figs. UX C ) an d (d), we show the phase diagrams for U12 = —0.7; the former is obtained 
by applying the GA directly to the original Hamiltonian of Eq. (Cp) while the latter is obtained 
from the effective Hamiltonian of Eq. ([8]). In the phase diagrams, there are only the MI 
phases with even total fillings, and the local state of the MI phase is |f , f )■ We see that 
the phase transition between SF and MI is entirely of the first order. Although the first- 
order SF-to-MI transitions for attractive inter-component interaction have been previously 



30] , our result that extends the 



found at commensurate fillings [l^] and in the hardcore limit 
region of the first-order transition is complimentary to the previous findings. In Fig. [l^c), it 
is shown that there are direct transitions between different MI phases even when t/U > 0. 
However, this is an artifact stemming from the fact that the approach fails to describe the 
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0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

zt/U zt/U 
(c) (d) 

FIG. 1: Phase diagrams on the (zt/U,/j,/U) plane for (a) Uu/U = 0.9, (b) Uu/U = 0.5, and 
(c), (d) U12/U = —0.7. The solid lines and the dashed lines denote the phase boundaries for the 
second-order and the first-order transitions. The phase boundaries in (a), (b), and (c) are obtained 
by solving Eq. ([3]) while the phase boundaries between PSF and MI in (d) are obtained by solving 
Eq. pop. The inset in (d) shows a magnified view that focuses on the PSF region between vacuum 
and the rit = 2 MI. In (a), At represents the width of the region of first-order transition in the 
i-axis and the dotted line is the trajectory used in Fig. [H 



PSF phase because of the lack of the second-order hopping process. Indeed, when the phase 
diagram is calculated from the effective Hamiltonian as shown in Fig.[T(d), the MI phases are 
separated from one another by thin but finite regions of PSF. Notice that although the MI 
regions in Fig. [H(d) are significantly larger than those in Fig. []Jc), this happens because the 
effective Hamiltonian is invalid for a relatively large value of zt/U > 0.1. Thus, combining 
the information of the two phase diagrams, PSF is located in the regions of small hopping 
zt/U < 0.1 sandwiched between MI regions. 
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0.15 0.155 0.16 0.165 0.17 0.175 0.18 0.165 0.17 0.175 0.18 0.185 0.19 

zt/U zt/U 

(a) (b) 

FIG. 2: Condensate density p s = \<&i\ 2 + \&2\ 2 with rit = 2. (a): The filled circles, the open circles, 
and the open squares denote data for U12/U = 0.9, 0.8, and 0.7. (b): The filled circles, the open 
circles, and the open squares denote data for U12/U = 0.6, 0.5, and 0.1. 



B. First-order phase transitions 

Let us now focus on the first-order phase transitions from SF to MI with even total 
fillings for repulsive inter-component interaction. To confirm the first-order transition at 



commensurate fil 



model in Ref. [l4|, we show in Fig. |2]the condensate density p s = |$i| 2 + | < &2| 2 versus zt/U 
across the transition point for n t = 2 and several values of U12/U. In Fig. [2](a), we plot 
p s for the values of U12/U at which transition is of the first order. The transition point is 
determined as the crossing point of the energy of the MI state and that of the SF state. 
At the transition point, the condensate density exhibits a jump, which we define as A, and 
the magnitude of the jump depends on Uyi/U . We find that the transition point (zt/U) tT 
decreases with increasing U12/U. Notice that although we show in Fig. |2]Ja) only the case of 
n t = 2, we have checked that the presence of the first-order SF-to-MI transitions is a general 
feature in even total fillings. In Fig. [2]^b), we show p s for the values oiU^/U at which the 
transition is of the second order. In this case, p s changes continuously across the transition 
point and the phase transition occurs at zt/U = (zt/U) t v = 0.172 regardless the value of 

u 12 /u. 

In order to determine the parameter regions where the transition is of the first order, 
we plot in Fig. |3]^a) the jump A as a function of Uyi/U at the top of the n t = 2 Mott 
lobe. A is zero for Uyi/U < 0.65 while it is finite, i.e. the transition is of the first order, 
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ings, which has been pointed out for the effective two-component J-current 





0.6 0.65 0.7 0.75 0.8 0.85 0.9 
U 12 /U 



0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 
U 12 /U 



(a) (b) 

FIG. 3: (a) Jump of the condensate density at the transition point A as a function of U\ijU at 
the top of the n t = 2 Mott lobe, (b) The width of the region of first-order transition in the t-axis 
At, which is indicated in Fig. Q^a) is shown in the unit of t c as a function of UyijU ■ 



for U12/U > 0.65. A is peaked at U12/U ~ 0.93, and it reaches zero at U12/U = 1. 
Recall that the phase separation occurs when U12/U > 1, and the SF-to-MI transitions in 
the phase separated gases are of the second order, because they are equivalent to those of 
one- component bosons. 

We next consider the transition between SF with incommensurate fillings and MI. In 
Fig. |3]^b), we plot the width of the first-order transition region in the t-axis At/t c as a 
function of Uu/U, where t c is the hopping for the transition at the top of n% — 2 MI lobe. 
We find that At/t c increases with increasing U12/U. In the limit of U12/U — > 1, At/t c reaches 
unity; this means that the entire phase boundary becomes of the first order. However, At/t c 
suddenly drops to zero at U12/U = 1 because the transition becomes of the second order. 

We plot the total filling factor n t as a function of fi/U along zt/U = 0.112 line for 
U12/U = 0.9 in Fig. HJ We see that the total density exhibits a jump at the transition 
points, which is similar to the condensate density jump discussed above. For zt/U = 0.112, 
it is found that the density jumps are 0.082 at [i/U ~ 1.07 and 0.02 at fi/U ~ 1.63. For a 
fixed value of zt/U the density jump at the small chemical potential is always larger than 
that at the large chemical potential. It is worth noting that measuring the density jump is 
a possible option to observe the first order transition in experiment. In the presence of a 
trapping potential, the local chemical potential spatially varies such that the transition from 
SF with commensurate filling to MI can occur. Within the local density approximation, the 
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FIG. 4: Total density n t as a function of (i/U along zt/U = 0.112 line for U12/U = 0.9. The 
trajectory is indicated by the dotted line in Fig. [1] (a). 



density jumps emerge at the edges of the n t = 2 MI plateau. Given the fact that first 
order transitions in a system of Fermi- Fermi mixtures with population imbalance have been 



observed by measuring density jumps 701 ]. it is expected that this way may work also for 



the first-order SF-to-MI transitions in Bose-Bose mixtures discussed here. 



IV. EXCITATION SPECTRA 

In Sec. mi we presented the formulation for calculating elementary excitations of a Bose- 
Bose mixture in a D-dimensional hypercubic optical lattice. Having obtained the ground- 
state phase diagrams in Sec. IIH] we use the formulation to reveal properties of excitation 
spectra in this section. Henceforth, we assume D = 2 and the momentum of excitations to 
be kx ~~ ky — — k. 



A. MI phase 



We first consider the MI phase with even total fillings and analytically calculate the 
excitation spectrum. For the MI phase of this type, the variational parameters are given 



by fni,r 



ni,n w ri2,n 



and the superfluid order parameters $ a are equal to zero. Hence, the 



linearized equations of motion (jl~8l) take 8x8 matrix form. Substituting these parameters 
into Eq. ffl~8]) and solving it, we obtain the excitation energy, 



(±) 



e(k) 2 - e(k)U (An + 2) + U 2 ± (e(k) - U(2n - 1) - 2U 12 + 2/x) 



(33) 
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where the plus (minus) sign corresponds to the particle (hole) excitation. The particle and 
hole excitations of component 1 are degenerate with those of component 2, and this means 
that there are four excitation branches in total. There are other solutions of Eq. ( 1181) that 
are independent of k with the energies, 

tujm u m 2 = -^m-iimi - 1) + ^-m 2 (m 2 - 1) + U 12 m 1 m 2 - ymi - fim 2 

- —nx(ni - 1) - —n 2 (n 2 - 1) - U X2 n x n 2 + ym + [m 2 , (34) 

where m a is an non-negative integer other than (m l5 m 2 ) = (ni,n 2 ± 1) or (ni ± 1,72,2)- 
Each branch corresponds to multi-particle or multi-hole excitation and its energy is always 
positive. 

We next consider the MI phase with odd total fillings. In this phase, the variational 
parameters are given by f$ m = ± (S ni)int+1 )/ 2 5 n2tint ^ 1)/2 + 5 nit{nt ^ 1)/2 5 n2t{nt+1)/2 ) . Inserting 
these parameters into Eq. ffl8|) . we obtain the excitation spectrum by diagonalizing the 
Bogoliubov equation. For the case of the MI phase with n t = 1, we obtain 

(U-fi- e(k) ± tvjj) [Gu =F tvjj) {{U 12 - y - e(k) ±hw){U -y- e(k) ± Hu) - e(k) 2 } 

-e(k) (U 12 - n ± hu) {U - n ± ftJ)] = 0. 

(35) 

In Fig. [5l we plot the excitation spectrum of the MI phase with n t = 1. There we see 
four branches in this MI phase, which correspond to one hole excitation and three particle 
excitations. In contrast to the MI phase, the excitations are not degenerate. There are other 
multi-particle excitations, which we do not explicitly show here. 



B. SF phase 

We now calculate the excitation spectra in the SF phase by numerically diagonalizing 
Eq. (j!8p . In Fig. E^a), we plot the excitation spectrum for n t = 1 and Ui 2 = 0. There we see 
two dispersive modes and two non-dispersive modes. While the gapful dispersive mode near 
the MI transition corresponds to oscillations of the amplitude of the order parameters, the 



gapless dispersive mode is a phase-fluctuation mode called the Bogoliubov spectrum 71 ]. 



Since the system consists of two independent and equivalent one-component Bose gases, the 



dispersive modes agree with those for a one-component Bose gas 
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431 ] and each branch is 
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FIG. 5: Excitation spectrum in the MI phase with n t = 1, Uu/U = 0.9, t/U = 0.01, and fi/U = 0.4. 
The lowest branch is a hole excitation and the upper branches are particle excitations. 

doubly degenerate. In contrast, the non-dispersive modes are particular to the Bose-Bose 
mixture. 

We plot the excitation spectra for the repulsive inter-component interaction in Figs. Efb)- 
(d). It is clearly seen that finite Uyi splits degenerate modes into in-phase and out-of-phase 
modes. Increasing U\2, the two gapless modes repel each other more strongly than the 
amplitude modes. In Fig. E(c), we find that a level repulsion occurs around ka/n = 0.4 
and the mode that was non-dispersive at TJ\i = acquires fc-dependence. When XJyijU = 
1, a non-dispersive mode appears again in Fig. [6]^d). Furthermore, for Uyi/U > 1, the 
gapless out-of-phase mode has an imaginary part around k = 0, which means the dynamical 
instability toward phase separation. 

In Fig. [71 we plot the excitation spectra for Uu < 0. We again see the repulsion of 
gapless modes and the shift of other modes compared with that of C/ 12 = 0. One of the main 
differences from the case of U\i > is that the gapless in-phase mode exhibits dynamical 
instability, which leads to collapse of the mixture, for strong attractive interaction. In our 
calculations, the mixture collapses at \ijU = —0.79 for t/U = 0.07 and n t = 1. 

C. Around the phase transition 

In Fig. [HI we plot the excitation spectra for n t = 1, Uyi/U = 0.9, and several values of 
t/U across the SF-to-MI transition point t c /U ~ 0.0406. Notice that when we calculate 
the excitation spectrum of the MI phase, we use the same chemical potential as that at 
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the transitional point. When t/U decreases from SF into MI, the following two remarkable 
changes happen in the excitation spectra. One is that the in-phase amplitude (gapful) 
mode approaches the in-phase Bogoliubov (gapless) mode and that the amplitude mode 
becomes gapless at the transition point so that its low-energy part coincides with that of the 
Bogoliubov mode. In the MI phase, the two branches split into gapful modes, corresponding 
to the particle and hole excitations. This behavior is analogous to the case of one-component 



Bose-Hubbard model 



39 



41 



43 



3M45| . The other is that the out-of-phase mode becomes non- 
dispersive and reaches zero at the transition point. It remains non-dispersive and zero in the 
MI phase. However, this property is an artifact of our approach in which the second-order 
hopping process is neglected [13j. In Sec. IIVD} we show from calculations on the basis of the 
effective Hamiltonian that this mode has phonon-like dispersion, reflecting the counterflow 
superfluidity. 
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FIG. 7: Excitation spectra of the SF phase with n t = 1 at t/U = 0.07, U\i = —0.4 (a) and —0.7 
(b). 




FIG. 8: Excitation spectra along the rit = 1 line at U12/U = 0.9. The hopping amplitude is (a) 
t/U = 0.05, (b) t/U = 0.045, (c) t/U = 0.042, (d) t/U = 0.041, (e) t/U = 0.040557 ~ t c /U, and 
(f) t/U = 0.039. 



In Fig. [91 we plot the excitation spectra for n t = 2, Uu/U = 0.5, and several values of 
t/U. When t/U decreases in the SF phase, two of the gapful modes approach and their gaps 
descend toward zero. Moreover, the two gapless modes also approach. At the transition 
point, the low-energy parts of these four modes coincide. This reflects the fact that both 
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FIG. 9: Excitation spectra along the nt = 2 line at UvijV = 0.5. The hopping amplitude is (a) 
t/U = 0.05, (b) t/U = 0.045, (c) t/U = 0.044, (d) t/U = 0.0435, (e) t/U = 0.04289 ~ t c /U, and 
(f) t/U = 0.041. 



of the in-phase and out-of-phase superfluids disappear at the same time at the transition 
point. In the Mott phase, there are a non-dispersive branch and two dispersive branches, 
which were discussed in Sec. IIV Al 

In Fig. [TUl we plot the excitation spectra for n t = 2, Uyi/U = 0.9, and several values 
of t/U across the SF-to-MI transition. Because of the first-order transition, the excitation 
spectrum changes discontinuously at the transition point as clearly shown in Figs. ITOT d) and 
(e). In contrast to the cases of the former two cases, any gapful modes do not coincide with 
the sound modes at the transition. In the MI phase, there is a non-dispersive mode expressed 
by Eq. with (mi, m-i) = (2,0) or (0,2) below the particle and hole modes. This mode 
becomes dispersive when the second-order hopping processes are taken into account. 



D. PSF and CFSF phase 

In Fig. HU we plot the excitation spectra in the PSF and CFSF phases with n t = 1. We 
set the parameters t/U = 0.03 and U\i/U = —0.7 for the PSF phase and t/U = 0.03 and 
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FIG. 10: Excitation spectra along the nt = 2 line at UyijU = 0.9. The hopping amplitude is 
(a) t/U = 0.05, (b) t/U = 0.045, (c) t/U = 0.04, (d) t/U = 0.0378, (e) t/U = 0.0377, and (f) 
t/U = 0.035. 



Uu/U = 0.9 for the CFSF phase. While we see a gapless mode and a gapful mode in the 
PSF phase, only a gapless mode is present in the CFSF phase. Gapful modes are absent 
in CFSF because the local state consists of only two Fock states |1, 0) and |0, 1). It is well 



known that the local states have to consist of t 
mode to appear in the excitation spectrum 39|, 



iree or more Fock states in order for a gapful 



4l| . If one calculates the excitation spectrum 



of the CFSF phase with n t > 3, one obtains gapful branches in addition to a gapless branch. 
The gapless mode in each phase exhibits linear dispersion, which reflects superfluidity of 
pairs or anti-pairs. As was discussed in Sec. IIV C\ if we calculate the excitation spectrum 
for the same parameters by solving the linearized equations of motion (Tl8|) with the matrix 
elements of Eqs. (fl~9l) and (120]) . the gapless modes are non-dispersive so that the Landau 
critical velocities are zero. This means that PSF and CFSF acquire superfluidity through 
the second-order hopping process. 
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FIG. 11: Excitation spectra of the (a) PSF phase with t/U = 0.03 and U 12 /U = 0.9, and (b) CFSF 
phase with t/U = 0.03 and U 12 /U = -0.7. 

V. DYNAMICAL STRUCTURE FACTORS 



Dynamical structure factors of one-component Bose gases in optical lattices have been 
experimentally measured by using the Bragg spectroscopy techniques, and this measurement 



has led to the observation of the Bogoliubov mode in the SF phase |48j. In this section, 
we calculate the dynamical structure factors of Bose-Bose mixtures for SF, PSF, and CFSF 
phases to show that these phases can be distinguished through the Bragg spectroscopy. 

In the Bragg spectroscopy, one exposes a Bose gas to an oscillating external field with 
momentum k and frequency u, which induces density fluctuations of the gas, and measures 
the response to the fluctuations. Here we consider the following three types of density 
fluctuation with the form of Eq. (|26|) : (i) a density fluctuation to one-component = 
F\,% — n-n, (ii) an in-phase density fluctuation = F- n = + Tin, and (iii) an out-of- 



phase density fluctuation G™ 



nn — Tin- Substituting these fluctuations into the 



linear response formulae Eqs. ( 1291) and ( 131]) . we calculate the dynamical structure factors 
Sp(k,u). These types of density fluctuation can be implemented with use of component- 
dependent laser beams jsl, Q . 

Let us first consider the SF phase. In Fig. [121 we plot the dynamical structure factors for 
n t = 1, t/U = 0.07, and U12/U = 0.9. The parameters are the same as those in Fig. El^c). 
As shown in Fig. [T2T a). the out-of-phase (lower) gapless mode strongly responds to F\^. 
We can also see that the in-phase (upper) gapless mode and the high-momentum part of a 
gapful mode respond weakly to i*\ j. It has been pointed out that the low-momentum part of 
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FIG. 12: (Color online) Dynamical structure factors Sp(k, uj) of the SF phase with n t = 1 for 
t/U = 0.07 and U12/U = 0.9. Shown are the responses to the density fluctuation of one-component 
Fi i (a), the in- phase density fluctuation F- n (b), and the out-of-phase density fluctuation F° (c). 



gapful modes does not respond to density fluctuations also in the case of the one-component 
Bose- Hubbard model 41|. Fig. [T27 b) shows the dynamical structure factors for F- n . For the 



fluctuation of this type, strongly responding modes are the in-phase (upper) gapless mode 
and the high-momentum part of the gapful mode that is the fifth excitation branch shown 
in Fig. El^c). However, the out-of-phase (lower) gapless mode and the other gapful modes 
do not respond to the in-phase fluctuation. Fig. IT2Tc) shows the dynamical structure factor 
for F° nt . In this case, the out-of-phase (lower) gapless mode responds strongly. The fourth 
excitation mode also respond but this response is very weak. 

In Fig. [T31 we plot the dynamical structure factors for n t = 1, t/U = 0.07 and Uu/U = 
—0.7. As shown in Fig. [TBTa). the in-phase (lower) gapless mode strongly responds to the 
one-component fluctuation and no other mode significantly respond. The modes responding 
to the in-phase or out-of-phase fluctuations are the in-phase (lower) or out-of-phase gapless 
mode and the high-momentum part of a gapful mode. 

We next calculate the dynamical structure factors for the PSF and CFSF phases. In 
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FIG. 13: (Color online) Dynamical structure factors Sp(k, oj) of the SF phase with nt = 1 for t/U = 
0.07 and U12/U = —0.7. Shown are the responses to the density fluctuation of one-component F\^ 
(a), the in-phase density fluctuation F- n (b), and the out-of-phase density fluctuation F° ut (c). 
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FIG. 14: (Color online) Dynamical structure factors Sp(k, uj) of (a) the PSF at U12/U = —0.7 and 
t/U = 0.03, and (b) the CFSF phase at U 12 /U = 0.9 and t/U = 0.03. 

Fig. [Ml we plot Sp(k,u) with respect to the one-component density fluctuation We 
see that the gapless mode and no gapful mode respond to the one-component fluctuation 
for each phase. In the PSF phase, gapless mode responds also to F- n , but does not respond 
to F° ut . On the other hand, in the CFSF phase, gapless mode responds also to F° ut and 
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does not respond to F- n . These results illuminate the essential property of the PSF (CFSF) 
state that the out-of-phase (in-phase) motion is forbidden while the in-phase (out-of-phase) 
motion exhibits superfluidity. 

Since each phase exhibits different responses to the density fluctuations as shown above, 
measurement of the dynamical structure factors through the Bragg spectroscopy can be 
used to identify the phases. A similar suggestion has been made in Ref. 32|, where dynamic 
responses of the SF, PSF, and CFSF states to several types of trap displacement have 
been investigated. We note that our results qualitatively explain the dynamics presented in 



Ref. 



32). 



VI. CONCLUSIONS 



We have investigated ground state properties and excitations of Bose-Bose mixtures in an 
optical lattice using the Gutzwiller approximation. We have obtained the ground-state phase 
diagrams in the (fi/U, zt/U) plane and found that the SF-to-MI phase transition can be of 
the first order in a wide range of \ijJJ for the strong inter-component interaction with even 
total fillings. We have calculated the excitation spectra for the several phases by solving 
the linearized equations of motion. In the SF phase, we have found a few gapful modes 
that are regarded as amplitude modes near the SF-to-MI transitions. When t/U decreases 
toward the second-order SF-to-MI transitions at n t = 1 or 2, one or two of the gapful 
modes descend and their low-energy part coincides with the gapless Bogoliubov modes at 
the transition point. We have also calculated the excitation spectra in the PSF and CFSF 
phases, and shown that there is the gapless mode that supports superfluidity of pairs and 
anti-pairs, respectively. Moreover, we have computed the dynamical structure factors by 
means of the linear response theory. We have studied the responses to three types of density 
fluctuation for each quantum phase. We have shown that branches which respond to the 
density fluctuations are different in each quantum phase. These results suggest that the 
quantum phases should be identified by measuring the responses to the density fluctuations. 

One of the most characteristic properties of first-order transitions is the hysteresis, where 
the phase transition point depends upon its history, namely from which phase we start 
changing parameters to cause the transition. Recent theoretical work has studied quantum 
phase transitions in the hardcore Bose-Hubbard model with long-range interactions on a 
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triangular lattice and predicted that when the chemical potential is varied, the hysteresis 
associated with the first order SF-to-MI transition exhibits an anomalous behavior in which 



a standard hysteresis loop is not formed 72J. It will be interesting to study the hysteresis in 
the two-component Bose-Hubbard model in order to reveal whether or not such an anomalous 
behavior is present also in the absence of long-range interactions. 
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Appendix A: Linear response theory within the Gutzwiller approximation 



For Simplicity, we use the contraction basis \n) = 1 711,712) and wavefunction = 

YliY^nf™ \ n )i * n this appendix. The equation of motion for the coefficient fn\r) 

can be derived by imposing the stationary conditions on the effective action S = 
d 



ih£-H- H pert (r) 



) with respect to f$*'- 

Ef® (XG e~ ilJJT e Trr + X*G* e iLdT e VT ) 



(Al) 



Here, we define E = (^>\H\$). H pert {r) and G rntn have been defined as ( |26|) and ( l27l) . 
respectively. In order to consider the small amplitude oscillations, we assume the following 
form: 



(A2) 



= [/« + 5/«(r)_ 

where fn is the stationary solution in the absence of the external perturbation, which 
satisfies 

dE 



Miff 



Oft 



{%)* 



(A3) 



A i )_f{i) 
Jn — Jn 



Thus, the linearized equation of motion for 5fn is given by 



in 



(0 



d 2 E 



Jn — Jn 



m 



d 2 E 



m 



Jn — Jn 



(A4) 
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We assume the fluctuation in the following form: 



SfP(r) = V^e^ - v* n e^e^, (A5) 

where the Bogoliubov amplitude u i>n and v iin are time-independent. Inserting Eq. (|A5[) to 
(TA"4l) . we obtain 



(fku + ihr] + hwi) u i>n = ^ ^ 



x dfgdft 

Jm ^i^n^ 



d 2 E 



u 



t (»)_ A') 

J n J n 



J n J n 



(A6) 



and 



(ho + - ^Dj) v i>n =y^y^ 



Ef (*)* \ 
Jm ^r^vn^n- 



d 2 E 



Jn — Jn 



u 



A*) _ ;(») 

J n — J n 



(A7) 



We next expand the fluctuation 8fn in terms of the normal-modes as 

^n=E(^- S >S*)' (A8) 
v 

<V=£(Mn- 5 >ln*), (A9) 
v 

where uf'^v^J are the solutions of the Bogoliubov Eq. (ITS]) , and 1/ is the mode 



index. We assume the Bogoliubov amplitude u\\ and vf^ are normalized by 



M|2 _ lyWp 



lation X). £ n 
obtain 



1. With this normalization, they satisfy the orthogonality re- 



(")*„>') _ .,(")* ,("0 



<? V)I a From Eqs. (lA"6j) . flA7l) . M and (lA"9j) . we 



£ [ft (w + i»/ - w„) u$4, - * (<" + ii? + w„) 

V 

^\h{uj + ir] -uj v ) v$A v -h(u + ir] + w„) uf^Bl 



Ef^XG 
Jm "i^-'n,mj 

m 

Jm 'n^m.n- 



From Eqs. ( 1A10I) and ( 1A11I) . we obtain the solutions for A v and B v as 



MG|0)i 



-A; 



(0|Gt|v)* 



ft CJ + 177 — UJ V ft ^— ' U + 17] — U t 

(0\G\v)i x 1^ <i/|G*|0); 



1 (0GW , 1 



ft ^-^ o; + i?7 + cj„ ft t— 1 to + 177 + w, 



(A10) 
(All) 

(A12) 
(A13) 
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where we have introduced the compact notations of the matrix element as 



m,n u i n I i 



(A14) 

<°I G H* = E E (/£ >G m,n«S - wgGn.m/i?) , (A15) 
n m 

"I^IO^EERX^ - = (0|G^>?, (A16) 



W&Wh = E E (/#^»«$* - v^Cr^fW-) = (u\G\0)l (A17) 

n m 

Having obtained the linearized solution for o/A , we now calculate the fluctuation of the 
physical quantity associated with Ff. 



where 



(Fi) = E f^* F n,mf^ « (£>0 + W>, 

V- 1 i/O — / j J n ± n,mj n i 

n,m 



Using Eqs. (1X5]) . flA8|) . (IA9|) . ( fAT2|) . and ( TAT3]) . we obtain 

JFI^^HGIO), (0|GH>|F|0) 



1/ 1' L. \ 



w + i?7 - 



W + 17] + UJ U 

X 3 e iujT e 



LO — 17] — UJ V UJ — 17] + UJ V 

This can be written in terms of the linear response function as 



where 



■-Y 



{Q\F\v) i {v\G\Q) j ^\G\v)^v\F\% 



H \ U + 171 — U) v U + 17] + bJ v 

Using the Fourier transformation, Eqs. (1A22|) and f ]A23j) are written as 

S(F t ) = [XF, G (k, oo)e-^ + XfM^ -fc^e*] e ikr *A k , 



(A18) 



(A19) 
(A20) 



(A21) 



(A22) 



(A23) 



(A24) 
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and 



XF,c( k > W ) 



1 



E 



(0|F|A)(A|G|0> (0|G|A>(A|F|0> 



(A25) 



h 



u + irj- U! KX uj + iri + w kiA 



A 
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